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CRITICAL THRESHOLDS IN EULER-POISSON EQUATIONS 


SHLOMO ENGELBERG, HAILIANG LIU, AND EITAN TADMOR 

DEDICATED WITH APPRECIATION TO CIPRIAN FOIAS AND ROGER TEMAM. 

Abstract. We present a preliminary study of a new phenomena associated with the 
Euler-Poisson equations — the so called critical threshold phenomena, where the answer 
to questions of global smoothness vs. finite time breakdown depends on whether the 
initial configuration crosses an intrinsic, 0(1) critical threshold. 

We investigate a class of Euler-Poisson equations, ranging from one-dimensional prob¬ 
lem with or without various forcing mechanisms to multi-dimensional isotropic models 
with geometrical symmetry. These models are shown to admit a critical threshold which 
is reminiscent of the conditional breakdown of waves on the beach; only waves above 
certain initial critical threshold experience finite-time breakdown, but otherwise they 
propagate smoothly. At the same time, the asymptotic long time behavior of the solu¬ 
tions remains the same, independent of crossing these initial thresholds. 

A case in point is the simple one-dimensional problem where the unforced inviscid 
Burgers’ solution always forms a shock discontinuity except for the non-generic case 
of increasing initial profile, u'q > 0. In contrast, we show that the corresponding one 
dimensional Euler-Poisson equation with zero background has globa l smo oth solutions as 
long as its initial (po,uo)- configuration satisfies Ug > —■\/2kpo - see ( ^.ll[ ) below, allowing 
a finite, critical negative velocity gradient. As is typical for such nonlinear convection 
problems one is led to a Ricatti equation which is balanced here by a forcing acting 
as a ’nonlinear resonance’, and which in turn is responsible for this critical threshold 
phenomena. 


AMS subject classification: Primary 35Q35; Secondary 35B30. 
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1. Introduction 


The Euler-Poisson equations 

(1.1) pt + V-{pu)=Q, a: G M”, t G M"*", 

(1.2) {pu)t + V • {pu ® u) = kpVcj) + viscosity + relaxation, 

(1.3) Acj) = p + background, x G 


are the usual statements of the conservation of mass, Newton’s second law, and the Poisson 
equation dehning, say, the electric held in terms of the charge. Here A; is a given physical 
constant, which signihes the property of the underlying forcing, repulsive if /c > 0 and 
attractive if /c < 0. The unknowns are the local density p = p{x,t), the velocity held 
u = u{x,t), and the potential 4> = 4>{x,t). 

This hyperbolic-elliptic coupled system 0)-(IID describes the dynamic behavior of 
many important physical hows including charge transport [^0| , plasma with collision |12], 
cosmological waves iQ and the expansion of the cold ions |llj| . Systems (0)- o also 
describe the evolution of a star regarded as an ideal gas with self-gravitation {k < 0) [T^ . 
The case of repulsive forces {k > 0) ||^ is relevant for plasma physics. These equations 
may be obtained from the Vlasov-Poisson- Boltzmann model by setting the mean free 
path to zero p. We would like to point out that the Euler-Poisson equation is closely 
related to the Schrodinger-Poisson equation via the semi- classical limit and the Vlasov- 
Poisson equation as well as the Wigner equation. Such relation has been the subject of a 
considerable number of papers in recent years; we refer to Q, |^] and references therein 
for further details. 

There is a considerable amount of literature available on the global behavior of Euler- 
Poisson and related problems, from local existence in the small Ll^-neighborhood of a 
steady state |]J, 1,1 to global existence of real solutions with geometrical symmetry , 

the relaxation limit for the weak entropy 
for isothermal case. 


for the two- carrier types in one dimension 
solution, consult 0 for isentropic case, and 


case, 

For the question of global behavior of strong solutions, however, the choice of the initial 
data and/or damping forces is decisive. The non-existence results in the case of attractive 
forces have been obtained by Makino-Perthame [^], and for repulsive forces by Perthame 
|^|. For the study on the singularity formation in the model with diffusion and relaxation, 
see [|^. In all these cases, the finite life span is due to a global condition of large enough 
initial (generalized) energy, staying outside a critical threshold ball. Using characteristic- 
based methods, Engelberg gave local conditions for the hnite-time loss of smoothness 
of solutions in Euler-Poisson equations. Global existence due to damping relaxation and 


with non-zero background can be found in 
relaxation the global existence was obtained by Guo 



For the model without damping 
assuming the flow is irrotational. 


His result applies to an iL 2 -small neighborhood of a constant state. Finally we mention the 
steady solution of non-isentropic Euler-Poisson model analyzed for a collisionless plasma 
in 1171 and for the hydro-dynamic semiconductor in |]^] — their approaches are based on 
the phase plane analysis. 

In this paper we present a preliminary study on a new phenomena associated with the 
Euler-Poisson equations — the so called critical threshold phenomena, where the answer 
to the question of global vs local existence depends on whether the initial configuration 
crosses an intrinsic, 0(1) critical threshold. Little or no attention has been paid to this 
remarkable phenomena, and our goal is to bridge the gap between previous studies on the 
behavior of solutions of the Euler-Poisson equations in the small and in the large. To this 
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end we focus our attention on the n- dimensional isotropic model, 
r'^nt + {nur'^)r = 0, r > 0, 


(1.4) 


p{ut + uur) = kpcpr + viscosity + relaxation, 
{r’^4>r)r = + background, u = n — 1. 


It is well known that finite time breakdown is a generic phenomena for nonlinear hyper¬ 
bolic convection equations, which is realized by the formation of shock discontinuities. In 
the context of Euler-Poisson equation, however, there is a delicate balance between the 
forcing mechanism (governed by Poisson equation), and the nonlinear focusing (governed 
by Newton’s second law), which supports a critical threshold phenomena. In this paper we 
show how the persistence of the global features of the solutions hinges on a delicate bal¬ 
ance between the nonlinear convection and the forcing mechanism dictated by the Poisson 


equation as well as other additional forcing mechanism on the right of (1.4). 

In particular, the persistence of the global features of solutions does not fall into any 
particular category (global smooth solution, finite time breakdown, etc.), but instead, 
these features depend on crossing a critical threshold associated with the initial configu¬ 
ration of underlying problems, very much like the conditional breakdown of waves on the 
beach; only waves above a certain critical threshold form crests and break down, otherwise 
they propagate smoothly. See for example [^] for ion-acoustic waves with such critical 
threshold phenomena. 

At present, no rigorous results exist on the question just raised concerning the critical 
threshold phenomena in equations of Euler-Poisson type. In this paper we provide a 
detailed account of critical threshold phenomena for a class of Euler-Poisson equations 
without pressure forcing. We use these systems to demonstrate the ubiquity of critical 
thresholds in the solutions of some of the equations of mathematical physics. At the same 
time, we show that the asymptotic long time behavior of the solutions remains the same, 
independent of whether the initial data has crossed the critical threshold or not. 

We note in passing that in this paper we restrict ourselves to the pressureless isotropic 
Euler-Poisson equations. The existence of the pressure allows for additional balance, and 
we hope to explore the critical threshold phenomena for the general, possibly non-isotropic 
model with additional pressure forces in the future work. 

A simple example of an equations with a critical threshold is the ID unforced Burgers’ 
equation, ut + uux = 0. This equation describes the movement of particles that are not 
being acted on by any forces. The variable u{x,t) represents the velocity of the particle 
located at position x at time t. The global existence is ensured if and only if Uq > 0. Thus, 
the Burgers’ solution forms a shock discontinuity unless its initial profile is monotonically 
non- decreasing. In this case the finite time breakdown of the Burgers’ solution is a 
generic phenomena. In contrast we show below that the corresponding ID Euler-Poisson 
equation with zero background has global smooth solutions if and only if Uq > —\/2kpo (see 
(2.II) below ), allowing a finite negative velocity gradient. This is the critical threshold 
phenomena we are referring to. As is typical for such nonlinear convection problems, one 
is led to a general Ricatti equation which is complemented by a particular inhomogeneous 
forcing dictated by the Poisson equation. It is the delicate balance of the latter, acting as 
a nonlinear resonance, which is responsible for this critical threshold phenomena. 

This paper is organized as follows. In Section 2, we consider the ID Euler-Poisson 
equations with zero background, and show that the solutions of the corresponding Cauchy 
problem blow up in finite time if and only if certain local “threshold” conditions on the 
initial data are met. In this case the density and the velocity gradient are shown to 
decay at some algebraic rates. For this simple model we utilize both the Eulerian and 
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the Lagrangian description of the flows to investigate the critical threshold phenomena. 
We also discuss how the solution behavior depends on the initial data as well as on the 
coupling parameter k. The behavior of the velocity gradient is outlined for various cases 
classified by the relative size of the initial data. 

In Section 3, we present the critical thresholds for the ID model with nonzero constant 
background and the effects of the damping relaxation. In this case the solution oscillates. 
The oscillatory behavior of the solution is induced by the presence of nonzero background, 
and the oscillations can be prevented only by strong damping relaxation. The zero limit 
of the background is shown to be a kind of “singular limit”, which coincides with (close 
to) the zero background case studied in Section 2. We also justify a rather remarkable 
phenomena, namely, that the non-zero background is able to balance both the nonlinear 
convective focusing effects and the attractive forces, to form a global smooth solution 
subject to a critical threshold in the attractive case A: < 0. Moreover in this case the 
density converges to the vacuum exponentially fast in time. 

In Section 4 we continue the discussion for the ID model with viscosity, which is similar 
(though not identical) to the one found in the Navier-Stokes equations. We obtain an 
upper threshold for the existence of global smooth solution and an lower threshold for the 
finite time breakdown, and consequently, these imply the existence of a critical threshold. 
If the initial data happens to be below the lower threshold, the solution must breakdown 
in finite time. We see that the presence of a self-induced electric field as well as additional 
viscosity are not necessarily enough to stop the formation of singularities. 

In Section 5 we consider the Euler-Poisson equations governing the v + \ dimensional 
isotropic ion expansion in the electrostatic fluid approximation for cold ions. We show 
that for each model with given integer u > 0, there exists a critical threshold condition. 
Several issues which are also clarihed in this section include: 

(1) Expansion rate of the flow path (consult the recent work Q of Dolbeault and Rein 
in this context). 

(2) The large time behavior of the velocity; 

(3) The decay rate of the density as well as the velocity gradient; 

(4) Sharp estimate of blow up time when the initial data exceed the critical threshold. 
More precisely we provide the explicit form of the critical threshold for the planar case 

and 4-dimensional case. For other cases (including the cylindrical and the spherical case) 
we conhrm the critical threshold phenomena by establishing both the upper threshold 
for the existence of the global smooth solution and the lower threshold for the finite time 
breakdown. A key step in the proof is to introduce the proper weighted electric held, which 
is shown to be constant along the particle path. This fact combined with the momentum 
equation gives the decoupled equation for the how map, and a nonlinear resonance is 
responsible for the critical threshold phenomena. 

2. Critigal thresholds-ID model with zero bagkground 
We consider the ID Euler-Poisson equation of the form 

(2.1) pt -h {pu)x = 0, 

(2.2) ut + uux = —k(f>x = kE, 

(2.3) Ex = -(t>xx = P- 

Here A: is a given physical constant, which signihes that the underlying forcing kE is 
repulsive if A: > 0 and attractive if A: < 0. The unknowns are the density of negatively 
charged matter p = p{x,t), the velocity held u = u{x,t), the electric held E = E{x,t), 
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and the potential (j) = 4>{x,t). Equation (E3) is a statement of the conservation of mass, 
equation (2.2) is a statement of Newton’s second law, and equation (2.3) defines the electric 
field in terms of the charge. 

To solve the problem on the half plane (x, t) £ M x R'*' we prescribe initial data as 


(2.4) p{x,0) = Po{x) >0, Poe C'^(R), 

(2.5) u{x,0) = uo(x), uq G C'^(R), 


and we show that the solutions of (p.lj)-( |2^ ) with the above initial data break down in 
finite-time if and only if certain local “threshold” conditions on the initial- data are met. 

Set d := Ux{x,t), then dx{2^) together with (^) yield by differentiation along the 
characteristics. 


d' + = kp, 

p' + pd = 0, ' := dt + udx- 


Multiply the first equation by p, the second equation by d and take the difference. This 
gives 

^ d Y pd' — dp' j 

W = 

and upon integration one gets 


d 

P 


P{t) with [3{t) := kt + u'o/po- 


The decoupled equations for d and p now read 

( 2 . 6 ) d' + d^ = ^^d, 


(2.7) p' + /3(t)Y = 0. 

^From these equations we can obtain the explicit solution formula for d = Ux and p, 
respectively. We want to point out that it is this time-dependent factor /3(t) balancing the 
nonlinear quadratic term that is responsible for the critical threshold phenomena. Indeed, 
with -gY = 0 one has the usual blow up, d{t) = associated with the unforced Ricatti 

equation d' + d"^ = 0. The presence of a forcing of similar strength on the right of ( |2.6|) , 
kt+u' /po to ‘nonlinear resonance’, which, as we shall see below, prevents blow up 

at least above a critical threshold, such that —^/2kpo < u'q < 0. To highlight this fact we 
present the following general lemma. 


Lemma 2.1. Consider the ODE 

wt = a{t)w + b{t)w‘^, w{t = 0) = Wo- 
It admits a global solution 


w{t) = -^- j -, B{t) := b{t)e^o <Ddr 


1 -wo/o B{T)d 
provided the initial data, wo, is preseribed so that 


Wo 


[ B(T)d 

Jo 


!r < 1 for all t > 0. 
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Proof. Set 

v{t) = w(t)e-fo . 
Substitution into the above ODE leads to 

vt = 

Note that vq = wq. Then the solution can be written explicitly as 

u = r(;o(l - / B{T)dT)~ , 

Jo 

from which the lemma immediately follows. 


□ 


As an immediate application of this lemma we check the conditions for do and po so 
that the global existence of Ux and p is ensured. Lemma 2.1 applies to the above equation 
for d, (^), with b{t) = —1, a{t) = and wq = do = Uq. It follows that if Uq > —\/2kpo 
then global regularity for d is ensured, for 

ft ft .r k j k 

1 — Wo B{T)dT = 1 + Uq / 6-'° ^drdr = 1 + UQt + —pot^ > 0. 

Jo Jo 2 

Similarly, for the p equation (^) we have a = 0, 6 = —fJ and wo = po, and hence 

ft ft 

1 - Wo B{T)dT = 1 + po / P{T)dT = 1 + UQt + -pot^ > 0, 

Jo Jo 2 

provided the initial velocity gradient u'q remains above the same critical threshold —^/2kpo■ 

We now turn to an alternative derivation of this critical threshold. This Lagrangian- 
like approach will prove to be useful for more general cases. We start by appealing to 
“physical” considerations. As the system being described has no external forces acting on 
it, momentum ought to be conserved. In fact, multiplying (2.2) by p, multiplying (2.1) by 
u, and adding the resulting equations, one finds the momentum equation in its standard 
conservative form 

{pu)t + {pu^)x = kpE = {kE^ / 2) X, p = Ex. 

It follows that the momentum is conserved, pu{-,t)dx = const, provided the boundary 
terms vanish, i.e., pu^ —> 0 as x —> Too, and in particular, E^(oo) - E2( —oo) = 0. It is 
reasonable to require that the total charge, E{oo), be finite, and since Ex = p > 0, this 
implies that E(oo) = —E{—oo), for otherwise, p = 0. Thus the electric field is given by 


Eix,t) = ^(^J p{^, t)df,- j p{^, t) d^ 


Equipped with this expression of E in terms of the density p, we can obtain the explicit 
solution along the characteristic curves, x{a,t), parameterized with respect to the initial 
positions, x{a, 0) = a, 


( 2 . 8 ) 


fL 

dt 


x{a,t) = u{x{a,t)), x{a,0) = a. 


The momentum equation (^) tells us that 

d 


(2.9) 


dt 


u{x{a,t),t) = kE{x{a,t),t), 
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and hence the electric field remains constant along x{a,t), 




2^(x(a,t)) •+ J Pt{C,t)dC- 

2u(x(a, t),t)p(x{a, t),t) 

px(a,t) 


' x{a,t) 


pt{i,t)di 


L 


x(a,t) 


ip{(,t)uiC,t))xd^ 


= 0 . 


Physically, the constancy of the electric field along characteristics is clear: as no charge 
can cross trajectories, the amount of charge to the right and to the left of any given 
trajectory is constant. And since the electric field along a trajectory is half of the difference 
of these numbers, the electric field on any trajectory must be constant as well. With 
E(q:, 0) =: Eo,u(a,0) =: uq, and p{a,0) =: po, we find from (^) that 

u{x{a, t)) = Uq + kEot. 

This together with (|2^) yield 

(2.10) x{a,t) = a + uot + kEot^/2. 

In the inviscid Burgers’ equation (corresponding to the case A: = 0) the straight charac¬ 
teristics must intersect in finite time leading to finite time breakdown. Here the straight 
characteristic curves are replaced by the characteristic parabolas, which explain the critical 
threshold phenomena. Indeed, since u = uq + kE^t and Eqq, = poj we conclude 


Ux{x{a, t),t) = {uq + kpot)/ 


dx{a, t) 


u'q + kpQt 


, duo{a) 

Mo := 


da 1 + u'^t + kp^t"^/2' ' da 

Integrating the /?—equation, which we rewrite as ■^p{x{a, t),t) = —UxP, we find that 


p{x{a,t),t) = 


Po 


r(Q:, t) := 1 -|- Mq! -|- kpot'^/2. 


r(Q;, t )' 

Clearly, the positivity of the so called “indicator” function, r(a, t) := , is a necessary 

and sufficient condition for existence of global smooth solution, \ux\,p < Const., and 
consequently, higher derivatives are bounded. Thus the solution remains smooth for all 
time if and only if —Uq < \j2kpo. Conversely, if there are points at which this condition 
is not fulfilled, then p and Ux will blow up at finite-time. 

To sum up, we state 


Theorem 2.2. The system of Euler-Poisson equations (2.1)-(2.S) admits a global smooth 
solution if and only if 


(2.11) 'Wo('^) > —'v/2^/Oo(«)) Va G 

In this case the solution of /i2.]\) - ( jg. j|) is given by 

Po .. _ K + ^Pot 


p{x{a,t),t) = 


Uxix{a,t),t) = 


T{a,t) 


r(Q;, t) := 1 -|- MqA -|- kpot^/2 


T{a,t)' 

so that p ^ t~‘^ and Ux ^ t~^ as long as po ^ 0. If condition (2.11) fails, then the solution 
breaks down at the finite time, tc, where r(Q;,tc) = 0. 

To gain further insight into the behavior of the solution, we now turn to discuss the 
dependence of the solution on the relative size of do = u'q and pQ as well as the parameter 

k. 
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To be specific we consider only the behavior ol d = Ux for the repulsive forces /c > 0, 
since the solution always breaks down in the attractive case /c < 0. We first look at the 
dependence on do = Uq. Figures 1-3 describe the three different scenarios for the evolution 
of Ux depending on the relative size of do and po. 

(1) do > 0 (see Figure 1). There are two such cases: 

(i) if do > y/kpo, then d is decreasing and satisfies 


0 < d < do, d ~ 2/t. 
(ii) if 0 < do < y/kpo, then we have 


0 ^ d ^ dmax, djxiax — 


kpo 


^j2kpo - d^’ 


where dmax denotes its local maximum taken at At this time d' = 0, therefore te 
satisfies 


k‘^pf 


0+2 


fg -j- kpodo T do kpo — 0, d^nax — d(tg ). 


(2) —\/2kpo < do < 0 (see Figure 2). There are also two cases: 
(i)if —y/kpo < do < 0, then d 
and then attains its maximum at 


(i)if —\/kpo < do < 0, then d starts to increase and becomes zero at to = 


y/2kpo^^Q - do ^ + 

tg -- - - > to- 

kpo 


In this case we have 


do ^ d < djnax — 


kpo 


^/2kpo - dl 

(ii) —y/2kpo < do < — V kpo ■ In this case we have 

dmin ^ d < draax: d ~ 2/t, 


, d ~ 2/t. 


where 


dmin — d(x, tg ) — 


-kpo 


V^2^ 


_ rp 

“o 


dmax — d(^X,t^ ) — 


kpo 


^/2kpo - d^ 


and 

I _ ±v^2A:po - dl - do 

^ ~ kpo 

(3) do < —\j2kpo (see Figure 3). In this case the solution must break down at time 
t = t~. The blow up time tc can be obtained via F(a, t^) = 0, that is, 

± _ -do ± Y^d^ - 2kpo 

“ kpo 
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Eigure 1. 0 < do, A: > 0 



Eigure 2. -y/2kpo < do < 0, A: > 0 



Eigure 3. do < -\/2kpQ, k > 0 


Next we look at the solution behavior as the parameter k changes. Figures 4-7 display 
such changes for various choices of k. Again there are three possible scenarios; 
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( 1 ) 

( 2 ) 

(3) 


k > dg/po (see Figure 4). There are two cases: 

(i) if<i„>0, then 0 < d < 

(ii) if4<o,the„rf„<^5te^. 


^ < k < dl /pq (see Figure 5) 

(i) if do > 0, then 0 < d < do! 

(ii) if do < 0, then - 7 ^^ < d < 


fcpo 

^j2kpQ-d'^ 


0 < A: < ^ (see Figure 6). 

(i) if do > 0, then d is decreasing and 0 < d < do. 

(ii) if do < 0, then d < do and d starts to decrease and becomes unbounded at 
time t~. 

In closing, let us note the remaining cases of A:; if A: = 0 then we have the decoupled 
Burgers’ equation; and if A: < 0, the solution always breaks down. In either case, 
there is no critical threshold which yields global smooth solution for an admissible 
set of initial data. 










Critical Threshold in Euler-Poisson Equations 


11 



3. Critical Thresholds- ID models with nonzero background 


3.1. The basic model with constant backgronnd. We now consider the system 

(3.1) pt + {pu)x = 0, 

(3.2) ut + uux = —k(f>x = kE, 

(3.3) Ex — (f^xx — P c, 

with constant “background” state c > 0. Here we require that f^(p(0 ~ c) = 0. 
The presence of a constant background c in the Poisson equation changes the “physical 
situation:” we are now working in a universe which has a fixed background charge density 
of —c. There is also an equal amount of movable charge, p(x). This is (an approxima¬ 
tion of) the situation, for example, inside a metal or a doped semiconductor. The fixed 
background charge of —c corresponds to the fixed positive charge of an element which has 
had an electron stripped from its outermost shell. The movable charge corresponds to the 
electrons that have been liberated from the atoms. 

Using the Lagrangian-like approach, we solve the system of equations (3.1)-(3.3). For 


this system, one cannot expect the total momentum of the negatively charged particles 
to be conserved-they are being acted on by an outside force. From (^^) we find that 
E(x) = J^^ipiC) ~ c)d^. As the net charge in our universe is zero, we expect that the 
electric field intensity to vanish at x = Too, and hence we require that p{Eoo,t) = c. 
Likewise we require the particles be at rest at far field, i.e., u{±oo,t) = 0. This says that 
far from the origin our system is “properly charge balanced” and at rest. 

As noted above, ( ^ ) says that = kE with ^ denoting the usual differentiation 
along the characteristics, ■^x{a,t) = u{x{a,t),t). Using E{x) = f^^(p(0 ~ ^)d^, and 
following the same basic steps as above, we find that ^E = —cu. Combining these two 
results, we arrive at 

(3.4) 

yielding, for fc > 0, 


d?u 


u{a,0) = uo, ut(a,0} = kEg, 


u{x{a, t),t) = uq cos{V^t) H— -j= sin 

V CrC 
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Here, parabolas (corresponding to c = 0) are replaced by a different geometry of char¬ 
acteristics, where ^x{a,t) = u{x{a,t),t) implies 


x(a, t) = ^ sin(\/cfct) -|-^(1 — cosiVT^t)) -|- a. 

Vck c 


Expressed in terms of the indicator function, T{a,t) := Xa{a,t), given by 

r(a, t) = 1 -|- — -(cos{'/c^t) — 1), 

Vck c 

we proceed as before to find that Ux = Uajxa is given by 

u'q{q) cos(\/cfct) + sin(\/cfct) 


(3.5) 


Uxix{a,t),t) = 


T{a,t) 


Note that Ux = Ei/L we find from (^) that p' = —pux = —pVt/T which in turn leads to 

Po(a) 


(3.6) 


p{x{a,t),t) = 


r(a, t) 

Clearly, there is a global smooth solution if and only if r(Q;, f) remains positive. For this 
to hold, we note that there exists r such that F can be rewritten as. 


Po{ot) , 

T(a,t) = — + ^ 

and hence, F(q:, f) > 0 if and only if 


ck 


+ 


poia) 


-1) sin(t-l-r). 


Upia)'^ ^ ( Po{a) 


ck 


+ 


-1 < 


po{a) 


This is equivalent to the condition that |ug(a)| < k{2po{a) — c) for all a G M. 

We can summarize the case with repulsive force A: > 0 

Theorem 3.1. Consider the system of Euler-Poisson equations B)-IS with constant 
background charge c and the repulsive force k > 0. Then it admits a global smooth solution 
if and only if 

(3.7) l^^o(«)l < '\/k{2pQ{a) — c), Va G M. 

In this case, the density oscillates around the nonzero background charge c, and the velocity 
gradient does not decay in time. If condition (3^) fails, however, the solution breaks down 
at finite time. 

Remark 3.1. The above threshold is sharp, as shown in the phase plane (F,Ft). Actually 
integration of (^) along characteristic curves x{a,t) yields 

(3.8) x” -|- ckx = kEQ{a) + cka, 

and differentiation with respect to a, combined with Eq = po — c, leads to 

F'' -I- cfcF = kpo. 

Its energy integral then becomes 

^F'^ + c/cF^ - kpoT = ^u'l - k{po - c), 
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with a trajectory which is an ellipse centered at (po/c, 0). Here, even if initially Uq > 0, 
then after some time E^ can still become negative, with the possible breakdown due to 
intersection with the E = 0 line. Condition is the precise condition which rules out 
such a scenario. Geometrically, as c tends to zero, the center of the above ellipse moves 
to the infinity and the closed curve splits into the parabola we met earlier with the zero 
background case. In this limiting case, once Uq > 0, the trajectory will always run away 
from the ‘singular’ axis E = 0, see Figures 7-8. 



7. 

c > 
0 



8 . 

c = 
0 


Remark 3.2. Note that as c ^ 0, the critical threshold condition {3R_) tends to |uQ(a)| < 
A/2fe/9o (a). This is almost the same condition we found when we considered 0-0, 
the condition given in ( 2.11 ). Here, however, the part of the condition that involves Uo(a) 
appears with an absolute value even in the limit as c ^ 0. The absolute value appears 
because of the oscillatory nature of the solutions of our equation. When the derivative of 
the initial condition is too negative for some large value of t, the equation loses smoothness. 
One finds that as c ^ 0, this time tends to infinity. The passage to the case c = 0 is, 
therefore, a kind of “singular limit.” 

We now turn to the case of an attractive force, k < 0. Elsing the same basic steps as 
before we get the same characteristic description, as the case k > 0. Integration of 
this equation subject to initial data (x,x')(a, 0) = {a,uo{a)) yields 


E, 


x{a,t) = a H-h -- 


En 


Uq 


— yj—ckt 


+ x(-^ + 


Uq 


\/—ck 


/—ckt 


y/—ck) 2 

Therefore, using Eq^ = Po(q:) — c, we find the indicator function E(q;, t) := given 

by 


(3.9) E(a,t) = ^ + i 
c 2 


Po 

c 


1 -^- 


u'q 




+ 


1 


1 - 


Po 


+ 


Uq 


^/—ck 


/—ckt 


y/—ck) 2 

We conclude the rather remarkable phenomena, namely that the non-zero background is 
able to balance both the nonlinear convective focusing effects as well as the attractive 
forces, to form a global smooth solution subject to a critical threshold. 

Theorem 3.2. Consider the system of Euler-Poisson equations 0 - ( |g. j| J with constant 
background charge c > 0 and subject to attractive force, k < 0. Then, it admits a global 
smooth solution if and only if 

po{a) 


(3.10) 


Ur 


(«)>- !- 


y/—ck, Va E 
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In this case, the density approaches the zero exponentially in time, and the velocity gradient 
remains bounded uniformly in time. If condition l \3.10i) fails, then the solution breaks down 
in finite time. 


Proof. As argued in the proof of Theorem 3.1, we have that 

p{x{a, t),t) = po{a)/T{a, t), Uxix{a, t),t) = Ttia, t)/r{a, t). 

It is necessary and sufficient to show that the indicator function r(a, f) > 0 for all f > 0 
if and only if ( 3.1C|) is met. The necessity is obvious since otherwise if (|3.10|) fails then 
the second parenthesis on the right of (^) is negative, and hence T{a,t) would become 
negative for t large. 

For the sufficiency, there are two cases: either u'q{oi) = — (1 — ^°^^'^ )\/—ck, in which case 


r(a,t) = ^+(l-Lh 
c V c 

and obviously, r{a,t) remains positive for all time t > 0. For the other case where 
Ug > — ^1 — \/—ck, the possible zeros t* of r(a,t) are determined by 


Po\ 


(3.11) 


/—ckt* 


— Po 

c ' 


2po _ 1 I ^0 
c ^ ^ -ck 




/—ck 


A simple check shows that the quantity on the right is less than 1, that is, either < 0 
or no such real t* exists. Therefore T{a,t) > 0 for all time t > 0. □ 


Remark 3.3. Let us, again, note that with the zero background model, the repulsive force, 
k > 0, is necessary for the global existence of the smooth solution. When the nonzero 
background is being taken into account, we could still have the global existence even when 
the force is attractive, k < 0. The balancing effect of k and c can be observed clearly from 
the above results. However, in both cases, we hnd that there is a local “ critical threshold” 
condition on the initial data such that the solution remains smooth for all time if and only 
if this condition is met. 


3.2. A constant background model with relaxation. We consider a further modifi¬ 
cation of our problem (^), (|3.3|) , where (|3.2|) is now augmented by a relaxation term 


(3.12) 


u 

ut -|- uux = —kcpx -= kE — 

e 


u 


e > 0. 


We still require that f^(p(0 “ c) = 0. 

We are now working in a universe which has a fixed background charge density of — c. 
There is also an equal amount of movable charge, p(x). This is (an approximation of) the 
situation inside a metal or a doped semiconductor. The term —u/e is a “friction term”, 
which, as we shall see, causes solutions to decay. 

As before, we cannot expect the total momentum of the negatively charged particles 
to be conserved; after all, they are being acted on by an outside force. As the net charge 
in our universe is zero, we expect that the electric field intensity at x = Too will be 
zero. Consequently we require that p(±oo,t) = c and tt(±oo,t) = 0 so that Eix) = 
f-oo(p(0 - c) d^ ^ 0 as X ^ oo. 

Using the same techniques that we have used previously, we find that as a consequence 
of (|3.12|) we have. 


—u = kE — u/e, 
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and proceeding as before, 


(Pu 1 d 

-—r = —cku - -u. 

dP e dt 


Let e > and define n = ^ck - This guarantees that the solution of the last 

ODE will consist of damped sinusoids and decaying exponentials, where we find, 

//NX =1 f /X kEo — uo/( 2 e) . , . 

u{x[a, t),t) = e I tig cos[iJ,t) H- 

V M 

We see that as long as the solution of our system is well dehned, u{x, t) must vanish as 
t oo, and since ^x{a,t) = u{x{a,t)), we conclude that the characteristic curves have 
the form 


-t 

e 26 


x{a,t) = a + Eq/c - —EQCOs{^t) + 


2ecuo — Eq 
2e^ 




As we have seen already, solutions will cease to exist when the indicator function r(Q;, t) = 
^x(a, t) vanishes, for 


d 1 

uPx{a,t),t) = —u{x{a,t),t )——- 
oa i (a, t) 


t oo. 


It is easy to see that this leads to a critical threshold condition, namely that the vanishing 
of the following indicator function 


(3.13) 


r(a, t) = — + -— ( (c - Po) cos{pt) + 
c c 


e 26 


2ecrtQ - po + c . 


2 ep 




To make precise the condition on the initial data to have the global solution, we need to 
verify the first time when the local minimum of r(Q;, t) = -^x{a, t) intersect the ‘singular’ 
r = 0 axis as shown in Figure 10. 



Figure 9. Phase plane graph of (r,rt) for c > 0 with damping 


Thus, if we denote this first minimal time as t*{a), then the global solution exists if 
and only if 


r(a,f) > 0. 
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In other words, the above expression of r(a!, t) implies that there is a global in time solution 
if and only if, 


J (c - p„r^ + 

y /X 2eiJ. 


and this is equivalent to the condition that 
c- Po 


Uq + 


2 ec 


< 




{e^ - 1) + —(2po - c), Va G 


Let us now compute t*. By definition, t* = is the smallest time such that 

rt{a,t'^)=0, > 0, n G N, 

where 0 < < • • • < t"" ^ oo. ^From ri(a, t") = 0 we have 

2ecu'„-po+c 1 / N 

/^x- kn -^(c-Po) 




2epuQ 


(c - po)p + f X u'q + 2ek{c - Po) ’ 


p = ck - 

4g2 


2e 2e/i 

Substituting this t” into the expression of Ttt we find 


^ t sinfut"') 

Fttia^n = ^ 


-U'n 


+ ^(c-po)) 


which shows that t” is a minimum with Tttict,tn) > 0 provided sgn(uQ) = —sgn(sin(/xt”')). 

Therefore t* is uniquely determined by 

(3.14) 

0 < t* < ^, if Uq < min{0,2efc(/?o — c)}. 


* 1 
t = -tg 
P 


-1 


2epuQ 


_u'q + 2ek{c - Po) 


if 2ek{po — c) < Uq < 0, 


2il ^ ^ F’ 

T < t* < if Uq > max{0, 2ek{po — c)}, 


< r < ^, if 0 < u'o < 2eA:(po - c). 

Moreover, we can now find the asymptotic behavior of p{x, t) as t —> oo. Indeed, following 
the same procedures as before, we hnd that p{x{a, t),t) = with r(a, t) given in ( ^.14 ), 
and it follows that p{x, t) ^ c exponentially fast, as one would expect. We summarize by 
stating that 

Theorem 3.3. Consider the system of Euler-Poisson equations (3^), ( S.li ), ^3.3i} with 


a constant background charge c > 0, a repulsive force k > 0 and weak relaxation e > . 

Let the critical time t* be defined in (3. W- Then if at all points a G M, 


(3.15) 


u'oia) + 


c- po{a) 


2ec 


< \ k- 


4ce2 


Po(«) 


(et*A - 1) + 2/90(a) - c, 


the solution of ( 3.P ), ^3.3^) is smooth for all time. In this case, u{x,t) —> 0 and 

p{x,t) —> c exponentially as t ^ oo. Otherwise, if condition ( d.ldj ) fails, then the solution 
of fi3.l3i) , ([3.3l) loses smoothness in a finite time. 

Remark 3.4. Note that as e ^ oo we recover the local condition ( 3.10| ) for the case without 
relaxation. 

Finally to complete our discussion of all choices oi k, c as well as e, we turn to consider 
the non- oscillatory case with strong relaxation e < In this case we have 


X =kE- 
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Note that E' = —cu = —ex', and hence 

E = Eo{a.) — c{x — a). 

Combining the above two equations one finds that 

x" H —X + ckx = k{EQ + ca), 


and this equation together with the initial data, {x, x'){a,0) = (a,uo(a)), leads to the 
characteristics of the form 


En 


Eo,,l 1 X _ , J_'| I ^ 


xia, t) = a + ^ + e-- [(-^)(- - —) - — + ^) + 


c 2 4e^ 2// 


c 2 4e/i 2/i 


Upon differentiation with respect to a, one finds the indicator function r(Q;,t) = Xa{ce,t), 

r(“>‘) = ? + '■* - T><5 - 4^> - + 4^) + ll'"' 

Based on this formula we claim the following. 


Theorem 3.4. Consider the system of Euler-Poisson equations l{3.3f ) with 

a constant background charge c > 0, a repulsive force, k > 0 and a strong relaxation term, 
e < . If at all points a G M 


(3.16) 


Po{a). 


UQ{a) > min <{ 0, -(1-( \l ^ - ck +— 


1 


1 


2e 


then the solution of (^. i| ), ( [3.71 ( \3-3i) is smooth for all time. In this case, u{x, t) —>■ 0 
and p{x, t) ^ c exponentially fast as t —> oo. 

Proof. Expressed in terms of A = a = ^(1 — and 6 = (1 — po/c)-^ + ^ with 
u = — c'fc, we have 

r(a, t) = 1 - 2a + (o - + (a + 

Note that r(a, 0) = 1 and r(a, oo) = 1 — 2a = ^ > 0. It suffices to show that r(a, t) > 0 
for all time t > 0 under the condition ( 3.16 ). First, if a + 6 > 0, then since a < 1/2 it 
follows that 

r = (1 - 2a) + (a + 11 + 

[ a + b } 

> (1 — 2a) + (a + 6)6*-^“^^* 

= (1 — 2a) + min{2a, a + b}e^‘'~^^'^ 

> min{l, 1 — 2a} > 0. 


Second, if a + 6 < 0, then one has to rule out a possibility of a local minimum achieved at 
positive time t = t* > 0. At such time we would have Tt{a,t*) = 0, which gives 


t* 


J_ / (u + A)(a -b) \ 
2u °^\{u-\){a + b)) ’ 


t*) = 2u{a + 6)(u - A)e('^“^)** > 0. 


Since (u — A)(a + 6) > 0, it follows that if 6 > ^a then (u + A)(a — 6)/(u — A)(a + 6) < 1, 
and hence that t* is negative. Put differently, t* < 0 and hence r(a,t) > 0 for all t > 0, 
if a + 6 < 0 and b > ^a. In summary of these two cases, a sufficient condition for the 
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global existence is 6 > min 
the definition of a and b. 


{^a, —a}, which is exactly the same as (3.16) when recalling 


□ 


4. Critigal Thresholds- ID model with viscosity 


We consider the solutions of a parabolic- hyperbolic version of (2.1)-(2.3)-the modified 
viscous Burgers-Poisson equations: 

(4.1) pt + {pu)x = 0, 

(4.2) 

(4.3) Ex = (pxx — Pi 

with repulsive force A; > 0 (if the matter is treated as charged particles). We show that 
despite the presence of a parabolic term on the right of (^), these equations can still lose 
smoothness in finite time when the critical threshold condition is crossed. 


ut + uux = -k(f>x + ( — = kE + \ — 


P 


P 


Note that the parabolic term on the right of 


is similar to the one, UxxjPi 


found in the Navier-Stokes equations, ut + uUx = Uxx!Pi with a difference term of — ^ . 

In this context, both types of viscosity regularization terms admit similar behavior. As 
/? —> 0 then Uxx/Pi —and hence all tend to infinity. More importantly, the 


forcing effect in (|4.2|) leads to u J, — oo, which in turn leads by (O), that p —> oo, and then 
Uxx/Pi — and hence all tend to zero. Thus viscous terms tend to zero when 

they are most needed and the solution blows up with p] oo. It is this vanishing viscosity 
mechanism which allows the possible blow up above the critical threshold in Euler-Poisson 
system (^)-(]T^). 

Let d = Ux- Differentiate (4.2) with respect to x to obtain 

(4.4) dt + udx = -d^ + kp + {d/p)xx- 
The difference, px(4.4) —dx (0) yields 

(4.5) pdt - dpt + updx - upxd = kp^ + p{d/p)xx- 


Our goal as before is to estimate the ratio (3 := d/p, thus ^x(^ 
parabolic equation 

Pt T uPx — k Pxx/ P- 


]) yields the following 


We would like to show a /3—bound due to a maximum principle of the form, 

(4.6) inf /9(a, 0) + kt < P{x, t) < sup /9(a, 0) -|- kt. 

^ a 

The difficulty here is that /o(x, t) must die at infinity if one wants the integral of d, that 
is if one wants u, to be finite at Too and one wants P to be bounded away from 0. If p 
dies at infinity, then viscosity has a coefficient that blows-up at infinity. We make use of 
the following theorem to show that even in such cases we still have a maximum principle. 


Theorem 4.1 (maximum principle). If ut + f{x,t)ux = a{x,t)uxx, u{x,t) < D{x)E{t) 
where D{x) is sub-linear in x and E{t) is exponential in time; u, ut, Ux, Uxx G C, \f{x, t)| < 
d; and 0 < a{x,t); a{x,t) G C, then the solution of the equation obeys the maximum 
principle: 

u{x, t) < sup rt(x, 0). 

X 
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To the best of our knowledge there seems no maximum principle of this form available 
in the literature. We present a self-contained proof of this theorem in the Appendix. 

To proceed we assume at most sub-linear growth of Uxjp at far field in space and 
exponential growth of Uxjp in time; i.e. there exists 7 > 0 such that for |x| >> 1 


(4.7) 


Ux 


<< Const.e'’'* 


x\. 


Equipped with the above maximum principle we establish the following 

Theorem 4.2. Consider the system (4-1)- with the repulsive force A: > 0. Assume 

that the smooth solution is sought such that Ux/p is sub-linear in the sense of satisfying 
(O- Then if for any a G M 


sup 

X 


u'oix) 

Po{x) 


< -i 


I 2k 

po(a)’ 


then (u, p) ceases to behave nicely in finite time. If 

inff 

^ Po{x) \j po{a)' 

then {u, p) remains smooth for all time. In this case as t —> 00 




dix,t) ~ -. 


Proof. Theorem 4.1 shows that (3 satisfies a maximum principle as long as u £ C^, |tt| is 
bounded, p £ C'^, and Ux/p is sub-linear in x and exponential in t. We assumed that these 
conditions all hold, and we make use of the maximum principle to bound f3 by ( [4.61 ). 
Because (O) says that the derivative of p along the curve x{a,t) is —pd, that is, 

P = -pd = -/3(^)p^ 

where P{t) = d/p satisfies the bounds stated in (|4.6|) . By Lemma 2.1 one has 

Po{a) 


p{x{a,t),t) = 


1 -h po{a) /q /3(r)dr 


This gives 


po{a) 


< p{x{a,t),t) < 


Po{a) 


po{a){kt‘^/2-I swpx fi{x,0)t)-11 \ 1 j-: j pQ{^a){kf^/2-\-mix f5{x,t))t)-\-1 

Finally, as we have bounds for fl = d/p we find that d must satisfy: 

po{a){kt + inix (3{x,0)) ^ u( r f\ f\ <■ Po{a){kt + supx P{x,0)) 

PQ{a){kt^/2 I-mix P{x,0)t)-11 ~ ’ ’ ~ po{a){kfi/2sup^ fl{x,0)t)1 

as long as both of the numerators remain positive. Suppose that A: = 1, the initial 
data po = uq = —2arctan(x). (Note that the total mass is finite.) We find that 

/3(x, 0) = —2. Our maximum principle argument says that either d —> —00 as t —> 2 — \/2 
or (n, p) ceases to behave well sometime before t = 2 — \f2. The “weakest” way in which 
(3 can cease to behave well is for /3 = Ux/p to grow faster than sub-linearly in space. As 
(3{x, 0) = —2, we find the supremum and infimum of (3 are equal. The above inequalities 
lead to equalities until the maximum principle no longer applies. We find that (as A: = 1); 

p{a) 


p{x{a,t),t) = 


p{a){t‘^/2-2t) + l 
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As long as p{x, t) is continuous in time this implies that up to and including the time 

< 2 — \/2 at which (3 stops being sub-linear in space p{x{a,t),t) is completely known. 
Moreover, for large a we know that /9(x(a,0),0) « 1/a^. (As by assumption \u{x,t)\ is 
bounded, we know that x{a, t) does not change quickly.) Thus if a is large enough we find 
that p{x{a,t),t) 7^ p{a) ^ Xjo?. This is sufficient to guarantee that p{x,t) does not decay 
faster than \jx^. In other words, if p{x,t) decays like 1/x^, it cannot decay any faster at 
a later time. If Ux/p grows non sub-linearly, then Ux will die at a rate that is slower than 
l/(ln(x)x)—that is u will tend to infinity. 

We note that if we take initial data in which po{a){kt‘^/2 + sup^,/9(a:, 0)t) + 1 is ever 
zero-if we ever cross this critical threshold-then the solution must either blow up in finite¬ 
time (because the solution has an upper bound that tends to —oo in finite- time) or the 
solution must cease to be nicely behaved in finite-time (so that the maximum principle 
ceases to apply). In either case, we find that once the critical threshold is crossed, the 
solution no longer behaves nicely for all time. 

The estimates above show that as long as /c > 0 and the denominators do not ever 
become zero, d —> 2/t and p —> 2/{kt^). We see that the solution of the parabolic- 
hyperbolic equation tends to behave like the solution of the hyperbolic one in the long 
term. Similarly, we find that if fe = 0 and the denominators remain positive for all time, 
then: 


t ~ ~ t 


where ci < 1 < C 2 . This is also the behavior that one would expect without the parabolic 
term (i.e. the inviscid Burgers equation). Hence the proof is complete. □ 


We see that the presence of a self-induced electric field and of viscosity is not necessarily 
enough to stop singularities from forming. This is not terribly surprising. The viscosity 
used tends to zero as p —> oo-it is smallest when it needs to be largest. This, of course, is 
a problem with the viscous term in the Navier-Stokes equations as well. Additionally, as 
we are dealing with equations in one dimension, the electric field that we are dealing with 
is essentially the field due to an infinite sheet of charge. The field due to such a sheet is a 
constant throughout space-forcing two such sheets to collide is easy. There is no reason 
that such a force should prevent the collision of sheets. What we have succeeded in doing 
is showing how the behavior of the solutions of our sets of equations reflects the behavior 
of the physical systems they represent. 


5. Critical Thresholds-Multi-D model with Geometrical Symmetry 

Let us consider the Euler-Poisson equations governing the u + 1 dimensional isotropic 
ion expansion in the electrostatic fluid approximation for cold ions 


(5.1) 

+ 

II 

(5.2) 

Ut + UUr = —k(f)r = 

(5.3) 

{r''(j)r)r = -nr'", 


subject to the initial data 

(5.4) (n,tt)(r,0) = (no,'«o)(?’), no{r) > 0. 

Here r > 0 denotes the distance from the origin, and /c > 0 is a known physical constant, 
which takes into account the general scaled version (|5.3D of Poisson equation {r'^(t)r)r = 
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—Airqur’^. The unknowns are the local particle density n = n{r, t), velocity u = u{r, t) and 
potential </> = (p{r,t). The geometrical factor u takes values 0,1,2, for planar, cylindrical, 
or spherical symmetry. 


Our interest is to find a critical threshold criterion for the Cauchy problem (5.1)-(5.4). 
Several issues that will also be clarihed along the way include: 

(1) Expansion rate of the flow path; 

(2) The large time behavior of the velocity; 

(3) The decay rate of the density as well as the velocity gradient; 

(4) Sharp estimate of blow up time when the initial data exceed the critical threshold. 

5.1. Analytic solution along the particle path. Without pressure forces, the particles 
move without collisions in a self-consistent electrostatic field. The system contains only 
particle-path characteristics. We shall trace the time dynamics along these characteristics. 
To this end we introduce the following weighted quantities 

e := Er'^ = —(j)rr'', p = nr^. 


With these definitions, (|5.l|) becomes 

(5.5) pt + {pu)r = 0, 


and the Poisson equation (5.2), now reads 

(5.6) Cr = p. 

To avoid having a singularity at the origin we require e(0, t) = 0, hence 

pr 

e = / p{C,t)di- 
Jo 

In view of (^), this implies that e satisfies the transport equation 

(5.7) et + UBr = 0. 

To solve this equation we define the “flow map” r{a,t) : 

dr{a, t) 


(5.8) 


dt 


= u{r{a,t),t), r{a,0) = a. 


Denoting differentiation along this characteristic curve by ^ := the mass equation (5.7) 


and the momentum equation (5.2) yield 

(5.9) e = 0, 

(5.10) r':=u{r,t). 

We now solve the above system subject to the initial data 

(r,e,tt)|i=o = (a,eo(a),uo(a)) 

with a G M'*' parameterizing the initial location, and with the weighted mass eo{a) = 
/J” noiO^'^dC 

By (|5.9|) , e remains constant along the isotropic characteristics, e = eo(a), and (CT) 
then yields 


(5.11) 


= r(a,0) = a, r {a,0) = uo{a). 


This relation shows that each fluid particle, starting from position a with initial velocity 
uo{a), is influenced by a central acceleration Once the solution of this equation for 
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r(a, t) is known, the other dependent variables u, p can be determined accordingly. We 
proceed to study the solution of (^.11). To this end, we introduce the ‘indicator’ function 


r{a,t) := e/o“'-(K«u),r)dr_ 

The geometrical interpretation of r(a,t) will be clear from the explicit solution of the 
Euler- Poisson system (|5 . ![) - ( ^ . 3D , given in the following lemma which will play an essential 
role in our discussion. 

Lemma 5.1. Consider the Euler-Poisson equations ( \5.!^) subject to the initial data 

(no, wo) G C^iR+) X Ci(M+). Let r(a, t) he the flow map defined in ^5. SJ , then 

_ dr{a,t) 

^ ^ da ■ 

Moreover, the solution of fjS.jj) - (5.4) is given by 

dr(a, t) 
dt ’ 
no(Q;)a^ 


(5.12) 

(5.13) 

(5.14) 


w(r, t) = 
n{r, t) = 

Ur(r,t) = 


r’^r(a, t) ’ 

rt{a,t) 


r(a, t) 

r{a,t) = u{r{a,t),t), r{a,0) = a, Va G 


Proof. Along the particle path one has 

fL 

dt 

Differentiating this equation with respect to a gives 

dfd \ d d 

lii V ^ '^r{r{,a, t),t)—r{a, t), -^r{a, 0) = 1. 

Hence -^r(a,t) = T{a,t) for any t G M+. 

The mass equation along the particle path r{a,t) becomes 

p' + pUr = 0, 

and integration in time leads to 

p{r,t) = po{a)/T{a,t). 

The formula for Ur follows from the definition of r(a,t). 

Using the expression of the solution, (^.12|) - (|5.14|) , we conclude 


□ 


Corollary 5.2. The smooth solution to the Euler- Poisson equations (5A)-(5^) blows up 
in finite time, t = tc, if and only if one of the following equivalent eonditions is met: 

(1) Jq" Ur{r{a,T),T)dT =-oo; 

(2) T{a,tc) = 0; 

(3) There exists an a G M such that ^{ct,tc) = 0. 


To ensure the existence of the global regular solution, therefore it suffices to start with 
the set of prescribed initial data for which ( recall ^(a,0) = 1), 

c)t 

— (a, t) > 0 Vt G M+. 
da 
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5.2. The isotropic flow map. Equipped with the above relations we are in a position 
to study the isotropic flow map r{a,t), and the zeros of ra{a,t) which characterize the 
formation of the singularity. 

We begin with the solution of the isotropic flow map r(a, t) governed by 
summarize its behavior in the following Lemma. 

Lemma 5.3. The solution of d-dimensional problem, r" = keo{a)r~'^, with initial data 
(r, r')(a, 0) = (a, uo(«) > 0) *5 as follows. (We classify the different cases by the value of 
u := d — 1. 

• u = 0; the flow map is given by 
(5.15) r(a, t) = a + uo(«)i + — 

and the velocity is 


5.11). We 


(5.16) u{r,f) = uo{a) + keo{a)t. 

• u = 1: the flow map is 


(5.17) 


/ / ^^0 2 I 

r(r, t) = a exp ( + uqt 


with the parameter r given implieitly by the temporal relation 

exp + uqt^ di- 

In this case we have 

(5.18) + ‘^OiUot + keof^ < r{a, t) < a + uot + Oftlnt). 
The eorresponding velocity is given by 

(5.19) u{r,t) = h{a,t), 

with h{a,t) ~ uo + 0{tlnt) determined implicitly by the identity 



• n = 2: r = r{a,t) is given implicitly by 

(5.20) = cosh-i(^ - 1) - cosh- 

R R 

with 

(5.21) Q{a) := 

In this case we have 




(5.22) [a^ + Sa^uot + ^keot^]^^^ < r{a, t) < a + ^Uq + 

the velocity is uniformly bounded in time, and 

(5.23) lim u{r,t) = Q{a). 

t—>O0 
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Proof. We start with u = 0. A straightforward integration of r" = A:eo(a) combined with 
the initial position a and the initial velocity UQ{a) gives the same formula for r that we 
met earlier, ( ^.10 ), in the 1- dimensional case. 

We then turn to the 2-dimensional case, u = 1, where the flow map equation ( 5.11 ) 
reads 

(5.24) r" = keo{a)r~^. 

Starting at the location a with velocity uq, then the energy integral is 

(5.25) 
implying 

(5.26) r = \Ju^ + 2keo ln(r/a), uq > 0. 

Let r(a, t) be a dimensionless parameter such that for constant a 


~ = ^eo(lnr - In a). 


, dt 
dr = — = 


dr 


= -J-d 


\J Uq + 2keoln{r/a) 


T r^J vlq + 2keoln{r/a) keo 

Then the parametric solution of ( ^.26 ) is given by ( ^.17 ), i.e., 

/ keo 

r = a exp I — 

Here, the parameter r is determined by t 


r{T = 0 ) = a. 


I So 2 I 
I - r + uqt 


dt = rdr = a exp 


keo _2 


+ uqt dr, 


which gives 


t = a 


r (keo 

i (- 


+ no? d( 


We want to show that r ~ 0{t) for large t, with the tight bound given by ( ^.18 ). To show 
the lower bound, we rewrite the equation ( ^.24 ) as 

{rr'y = keo + > keo, 

and integration twice gives the bound on the left of ( 5.18| ). Combining this lower bound 
with r” = keor~^ yields 

r" < keo[oy + 2auot + 

and integration twice gives the upper bound shown on the right of (5.18). 

Finally we turn to the 3-dimensional case u = 2, where the flow map equation reads 

(5.27) r" = A:eo(a)r“2. 

Starting at location a with velocity uq, the energy integral is 

(5.28) 

which gives 




r = Q{a){l — Rfr)^^'^ for uq > 0, 


with Q and R given in ( 5.21 ). Integration in time once gives the implicit formula of the 
flow map ( 5.22 ), where the dependence of r on a will play essential roles (i.e., Tq, = 0) 
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in our later analysis. We conclude with the large time estimate ( ]5.22| ). The upper-bound 
follows from the fact r' < Q{a). /,From the flow map equation ( ^.27| ) follows that 

{r^/Sy = 2r{r'f‘ + keo > keo, 

and integration twice yields the bound on the left of (|5.22D. □ 


Remark 5.1. The above formula reveals quite different geometry of the isotropic flow path. 
In the case u > 1, the energy integral ensures the lower positive bound for r. In fact 
recalling ( 5.25|) and (^.28 ), 


(5.29) 

(5.30) 
one finds 

(5.31) 

(5.32) 


(r) = Uq + 2keQlnr / a, u = 1, 

(ry = uS + ^(a>-'-r'-‘'), 


r{a,t) > ae , i/= 1 
V- 1 


u > 1, 


r{a,t) > 


2keQ 




1 

' v-l 


V > 1. 


Moreover following the same procedure as for u = 2 we obtain the expansion rate of the 
flow path for general situation u > 2, 


(5.33) + (u + 1)q^UQ t + ^ ^ ^ < r{a,t) <a + a . 

2 V V “ 


t. 


In the ID case, u = 0, we saw that if the initial velocity is negative then particles can reach 
the r = 0 line in a finite time due to the quadratic form (|5.15 ). In contrast, for u > 1, 
these positive lower bounds imply that if the velocity is initially negative, the particle path 
r may decrease only for a finite time and then increase due to the positive acceleration. 
To avoid the technical discussions here and in what follows we restrict ourselves to the 
case rto > 0. 


Remark 5.2. The above results show that for the case u = 0,1 the velocity grows linearly. 
But for the case u = 2, the velocity is uniformly bounded and converges to a positive 
constant as time becomes large. 

To study the critical threshold phenomena one may utilize the ‘indicator’ function 
r = Tq,. This study is carried out in next section. 

5.3. Critical thresholds. Using the above flow map we provide precise conditions on 
the initial data such that either the solution remains globally smooth or it breaks down 
in a finite time. 

We start by revisiting the 1-dimensional case. 

Theorem 5.4 (Global existence of smooth solutions for planar case (u = 0)). The smooth 
solutions of 0-0 with u = 0 exist if and only if 


(5.34) 


rto(a) > — y^2fcno(a), Va G 
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In this case the solution is given by 

n(r(a, t),t) 

Ur{r{a, t),t) 


noja) _ 

1 + u^t + |no(a)i2 ’ 
Uq + knot 
1 + u'^t + |no(a)t2 


Proof. Differentiating the flow map equation r = keo with respect to a. we find that 

r" = kno{a), 

for eoa = Po = no- The definition of E gives (r',r)(t = 0) = (ttg, 1). Thus the correspond¬ 
ing energy integral is 


[r']^ — 2knoT = [u'o]^ — 2kno- 

The geometry of the trajectory implies that to ensures the positivity of T the initial data 
should satisfy either Uq > 0 or [uq]^ — 2A:no < 0 for the case u'q < 0, which yields (|5.34). □ 


The above formula immediately yields the following 


Corollary 5.5 (Breakdown of smooth solutions for planar case (u = 0)). The smooth so¬ 
lution to the Euler-Poisson equations (0-0 blows up in finite time if and only if the 
condition, u'fia) > —■\/2kno{a), fails, i.e., if 

3a e M“'“ s.t. Uo(a) < —^/2kno{a). 

In this case, the density n{r,t) and Ur{r,t) become infinite as t ] T, where the blow-up 
time t = tc, is given explicitly by 

tc ■-= 2/ sup i -Uq + \J{uq^ - 2kno 


Remark 5.3. Consider the ID equation with nonzero background and additional relaxation 
term, where the equation (5.2) is replaced by 

ut 3- uur = kE — ufe with E^ = {n — c). 


Hence the equation for r reads 

r" = kE -. 

e 

Note that E' = —cu = —cr', thereby E = Eo{a) — c{r — a). These equations lead to an 
‘indicator’ function T = satisfying 

r' 

3-h kcT = kpo- 

e 

Using the phase plane analysis one recovers, for the ID half space problem, a G the 
same results obtained for the 1-D Cauchy problem, a G M, consult Section 2-3. 


The rest of this section is devoted to the multi- D case, u > 1, where we confirm the 
remarkable presistence of critical threshold phenomena in the multidimensional problem. 
But precisely determining sharp critical thresholds is far from trivial. Even in the isotropic 
case, the implicit solution formula makes the final conditions on initial data rather cum¬ 
bersome. Thus, for example, differentiating ( 5.11 ) with respect to a yields 

r" = kpor-^ - A:ueor-(^+Tr, 


which is coupled with the flow map equation r" = keor It is difficult to find an explicit 
sharp threshold for the initial data that distinguishes between cases for which T remains 
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nonzero and cases for which it does not. Though some further tedious calculations may 
enable us to obtain a complex criterion for cylindrical case (u = 1) as well as the spherical 
case (u = 2), we do not perform these calculations. Instead we give sufficient conditions 
for upper thresholds on the initial data for the existence of global smooth solution, as well 
as the lower thresholds for the finite time breakdown. These confirm the existence of an 
intermediate critical threshold —which is the focus of our interest in this work. 


We start with the 2D case. 

Theorem 5.6 (Global existence of smooth solutions for the cylindrical case u = 1). A global 


smooth solution of Euler- Poisson equations with u = 1 exists provided the ini¬ 
tial data (uojno) with Eq = a~^ Jq satisfy 


k 

Uq > - [anoh{a) — Eq], Va G 

Uq 


(5.35) 

where h{a) is determined by 

(5.36) ka^nouo / ^ , n; f =1, Va G M+. 

Jo [uq + ^kEoarjY/^ 

Proof. Recall the energy integral 


r = \JUq-\- 2keoln{r/a), 


from which it follows that 


/ 


r{a,t) 


df 


= t. 


\/ui + 2keoln{f/a) 
Differentiating the above equality with respect to a, one has 

r(aT) .. 

= ^(a,t), 


where 


A(a, t) := 


^/u'o 2keoln{r/a) 

uqu'o + kpolnf/a — keoa~^ 1 

la [uq + 2keolnf,/a\^/'^ uq 


Noting that eo = aEo, we introduce M := [kEo — uoUo]/{kano) and rewrite A in terms of 
r] := ln{f,/a), as 


, ,2 rj-M 

A[a,t) = kano . 2 , n, f -1772 

Jo [ug + 2kEoar]Y/^ 


e^dr] + —. 
Uo 


We shall show that A, and therefore E, remain positive for all f > 0 provided ( 5.35 ) holds. 
W consider the case M > 0, since the complementary case M < 0 is trivial. Note that 
eoo = Po = noa. A simple computation involving (|5.26 ) gives 

dA ln{r/a) — M 

— = kauQ— - 


dt 


Uq + 2kEoaln(r/a) 


In view of the monotonicity of r in time (^ > 0) we see that A may achieve its unique 
minimum at t*, where At{t*) = 0, i.e., r{a,t*) = ae^, and this minimum is positive 
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provided ( |5.35 ) holds, for 

i. 2 r 

A[a, t ) = ka riQ 

Jo 


> ka^riQ 


I 


0 [ttg + 2kEoari]^/‘^ ^ uq 

n-h{a) 1 ^ 

- -^e^dr] H-= 0, 

0 [mq + 2A:Eoar/]3/2 uq 


if M < h{a), which is equivalent to ( 5.351) . Therefore the indicator function r(a, t) = 
A{a, t)-\/u^ + 2keoln{r/a) remains positive because A{a,t) > A{a,t*) > 0 for all t > 
0 . □ 


Condition ( 5.35 ) could be viewed as an upper threshold in the sense of providing a 
sufficient condition leading to global smooth solutions, though the permissible class of the 
initial data for global smooth solutions is clearly larger. However, the existence of the 
critical threshold can be ensured by combining this upper threshold with the following 
lower threshold for the finite time breakdown. 

Theorem 5.7 (Breakdown of smooth solutions for cylindrical case u = 1). The smooth so- 

with u = 1 breaks down in finite time if 


lution to the Euler- Poisson equations (5.1)- 
the condition, Uo(a) > — y2fcno(a), fails, i.e., if 

(5.37) 3q; G M’*' s.t. UQ{a) < —\/2kno{a). 

Proof. Using the parametric form of the flow map given in ( ^.171) we evaluate the ‘indicator’ 
function T via 

r(Q;,t) = - rrjtrta. 

/,From (5.17) we see that 


Tq, = ( 1 + a;(|poT^ + Uqt) ) exp 


kcQ 2 , 

— T +UoT 


r^ = a [keQT + uq) exp 


keo 2 I 
—r +uor 


ia = ^ exp 


—f, +uof, 


di, 


keo 2 I 1 
—r + uqt ) = r. 


k 


tr = aexp 

Expressed in terms of 
(5.38) 

and r = r(r) = aexp[^r^ + uqt], the ‘indicator’ function can be rewritten as 

Jo^ rjOHQdQ 


h{T) := 1 + au'^T + -apor'^ 


V = 


a 


6 (r) - (no + keor)- 


r{T) j ■ 


Note that 6(0) = 1. The quadratic form of (5.38) implies that there must exist a parameter 
T* such that 6(r*) = 0 provided (|5.37|) holds. At this time T becomes negative because the 
nonlocal term r(^)b(^)d^ stays positive. This combined with the fact r(0) = 1 ensures 
that there must be a finite time t = t* such that T{a,t*) < 0. Hence T must vanish at 
finite time t = tc < t*. This completes the proof. □ 
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We conclude with the 3-dimensional case, stating the lower threshold for finite time 
breakdown. 

Theorem 5.8 (Breakdown of smooth solutions for spherical case u = 2). The solution of 
Euler-Poisson equations for u = 2 blows up in finite time if the condition, 

Uq > -:^[ano - Eq], fails, i.e., 


(5.39) 


k 

3a E M"*" s.t. Uq < -[ano — Eq], 

Uq 


Proof. Instead of using the implicit formula of the flow map r in ( 5.20 ) we introduce a 
dimensionless parameter r such that r may be rewritten in terms of the parameter r E M"*", 

R 

r = —[1 + cosh(r + ro(a))], 

R 

i + sinh(r + ro(a)) - sinh(ro(a))], 

where Q and R are given in (^.21|), and tq is determined uniquely by 


cosh(r„(a)) = - - 1 = —, 
Thus the ‘indicator’ function E is determined by 

r(a,t) = - rr/trta, 

with 

R R 

i’q = -^[1 + cosh(r + ro(a))] —sinh(T + To{a))toa, 

R 

rr = —sinh(r -k To(a)), 

R 

= ^[1 + cosh(r -k To (a))], 


tr = 


R \ R 

—j [t + sinh(r + ro(a)) - sinh(ro(a))] -k ^[cosh(r + ro(a)) - cosh(ro(a))]roa. 


Expressed in terms of 


A{t) := 


sinh(r + tq) 


1 + cosh(r -k To) ’ 


we have 


r(a, t) = 


1 -k cosh(r + To) 


-QAir) (R/Ql 


^Ra + RA{T)Toa 
T — sinh(ro) 


+ ' 4 ( 01 + 5 ( 1 - 


1 + COsh(To) 




1 -k cosh(r + Tq) j Q \ 1 + cosh(r + tq) , 

To show the breakdown of solutions it suffices to show that as r becomes large E becomes 
negative since E(a,0) = 1 > 0. Note that as r ^ oo one has A(t) —> 1, and therefore 
the limit of the sum in the bracket {• • • } combined with eo = and Cq = po = 

becomes 

R R 

gQa = + Ho^no - Eo)] . 
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Hence the ‘indicator’ function L would become negative for large time whenever the con¬ 
dition ( p.39D is satisfied. Therefore there must be a finite Tc, also a finite time tc, such 
that T(a,tc) vanishes. At this time, the solution breaks down. □ 

Remark 5.4. The above lower threshold for breakdown, u = 2, can also be verified in an 
alternative way as a particular case of the more general situation u > 2 . Indeed, the 
new feature for v > 2, is that the velocity tends to a constant for large time t’s, which is 
evident from the energy integral (^.30D, 


'Z I l\2 z I 
u ={r) =Uq + 


2keo 




l-u\ 


while noting that, say by 
constant value 




iz-1 

0 as t ^ oo. Hence, the velocity approaches the 


u{t) ~ (5(a), (5(a) = ^^0 + 

which in turn implies that 

r ~ Q{a)t. 

Consequently, if the following critical condition fails. 


2keQ 

v-l 


-a 


l-v 


Uq> - 

Uo 


noa 

iz-l 


— En 


i.e., if there exists an a G M'*' such that Qa < 0, or explicitly, that 


(5.40) 


k 

3a G s.t. u'n < - 

Uo 


nga 

u-1 


-En 


then two particle paths must collide at large time. Observe that this critical condition for 
u = 2 coincides with (5.3£). 


We conclude with the general v >2 case, discussing the upper threshold for the existence 
of global smooth solution. 

Theorem 5.9 (Global existence of smooth solutions for general v >2 cases). A global smooth 
solution of the Euler-Poisson equations (5.1)-( ^^) with u > 2 exists provided the initial 
data (uo, no) with Eo = a~'^ no{f,)f,'^d^ is prescribed such that for all a G K"'' 


(5.41) 


Uo > -- 


A 

Uo 


noa'^hp 


u-l 


Here, hp is determined by 

^ 2 ) kuonoa hp{a) - tj 


L 


— En 


{a^ — rf) ^-''drj = 1. 


(U - 1)2 [^2 + ^r?]3/2 


Remark 5.5. Since, as we shall see below, hp < a^ we conclude that the lower threshold 
(5.40) is indeed smaller than the upper threshold in (5.41). 


Proof. Recalling the energy identity (|5.30D, 


/ /\2 z z I 
(r) =u =Uo + 


2keo 

u-1 


(a^-^^-E-^) 


l-u\ 


dr 




dt. 


we have for tto > 0 
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Integration yields 


f‘r(a,t) 


di 


= t. 


Un 




I-iy', 


Differentiating the above equality with respect to a leads to 

r(a, t) 


Ufi + 


2A:eo 

v—1 


a 


l-u 


— r 


l-v 


= B{r,t) 


where 


B{a,t) := [ 
J a 


U0U0 + i^PoOi^ "-keoa 1 


which can be rewritten in terms of := (u — 1 ) ^^° and r] := ^ ^ 


B{a, t) = 


kuQa 


f 


kriQQ.^ 

y- Ml, 

[,„2 + 


as 


1 ^ 1 

[a — y) ^-•'dy H-. 

uo 


In remains to show the positivity of B for all t > 0 provided (^.41|) holds. If Mi, < 0, then 
it is easy to see B > 0 for all t > 0. We now consider the case Mi, > 0 by checking the 
positivity of the possible minimum of B. Note that eoa = Po = noa’^ and eo = Eqq^. A 
straightforward calculation involving (5.3C) gives 


dt 


knoa 


- M„ 




/,From the monotonicity of the flow map ^ > 0, it follows that there exists a time t* such 
that 

^l,=,.= 0 , ^(‘-«•)>() for 

and at this time r{a,t*y~'' = — M„. Therefore we have 

knoa y — Mi, 


B{a, t) > B{a, t*) = 
knoa I 


(U- 1)2 70 + 


y — h„(a) , 1 ._!±_ , 1 

>7 -/ —-- TTT - —y)^-"dy-\ -= 0 , 

{v-iyJo [u2 + f^7?]3/2 ^ ' UO ' 


provided for all a G M"'', M„{a) < hy{a) < a 


l-u 


i.e 


, (^.41 ), with E" defined in (5.42). 


Hence r(Q;,t) remains positive for all t > 0 once the initial data remain above the upper 
threshold ( 5.41 ). □ 

The above upper and lower thresholds for the cases u > 1 confirm the existence of 
an intermediate critical threshold, through we do not provide the explicit form of the 
critical threshold. One case in which we can precisely compute the critical threshold is 
the 4-dimensional (u = 3) isotropic case, which is given in the following 

Theorem 5.10 (Critical threshold for 4-dimensional model (u = 3)). The global smooth 

with u = 3 exists if and only if 


solution of Euler-Poisson equations (5.l)-( 

{auQ — uo)^ < 4a(^^^^ — Eq), Va > 0 and Mo(0) > 0. 


(5.43) 


2 
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In this case the velocity is given by 

auQ + [uq + kaEojt 


(5.44) 


u{r,t) = 


yj^ 2auot + [mq + kaEo]t‘^ 
and the density is given by 
(5.45) 


diiQ + taEo 3s t T oo, 


[a^ + 2auQt + (uq + kEoa)E][a + (uq + augjt + (uqUq — /cEq + |noa)i^] 

Proof. As argued before it is sufficient and necessary to show that the threshold condition 
(5.43) ensures the positivity of the indicator function E for all t > 0. Let us first solve the 
flow map equation (|5.11|) with u = 3, i.e., 

r" = keor~^, 


r(0) = a, r'(0) = uq. 


Its energy integral is 


[r^ = Uq + ke^a ^ — ke^r 


-2 


-2 


= ui + ke^a ^ — rr”, 


where r" = ke^r ^ has been used in the last equalities. Rewriting this relation leads to 


= [r'f + rr” = Uq + ke^a 


and integration twice gives 


Hence 

(5.46) 


+ 2auot + [ug + ke^a ‘^]t^. 
r = yj + 2auot + [uq + keoa~‘^]t‘^, 


which as u = -^ and eo = Ega'* gives the velocity (15.441) . The corresponding ‘indicator’ 
function is 

rlrv t) - — - ^ [<^uoyt + + kaEoft'^ 

dot + 2auot + [ug + /caEgjt^ 

_ a + [ug + au'olt + [uqUq - kEo + 

\JoP P 2auot + [u'q + fcaEg]f2 

Note that at the origin a = 0, r(0,f) = 1 + Ug(0)L These explicit formula imply that 
r(a,f) remains positive for all f > 0 once 

[uq + aug]^ < 4a(MgMg — kEo H- 2 ^^’ ^ ^ ^o(O) — 

This is equivalent to ( ^.43 ). The solution ( 5.45|) fol lows from the above explicit expression 
of r(Q;,f) when recalling the general formula (|5.13| ). □ 

The above critical threshold result enables us to claim the following 

Corollary 5.11 (Breakdown of smooth solutions for the case u = 3). A solution of the Euler- 
Poisson equations (B-(0 with u = 3 blows up in finite time if and only if condition 
{ 5 . 4 ^ fails, i.e., 


3a e M^s.L (oug — ug)^ > 4A:a( 


ngo: 


— Eg 


2 




























Critical Threshold in Euler-Poisson Equations 


33 


In this case, n{r,t) and Ur{r,t) become infinite as t ] T, where the blow-up time is given 
explicitly by 


tc := 2/ sup I -Mg - — + - 4 / (aug - - Eb) 1 . 

{ a a \ 2 j 


6. Appendix 

Proof of the maximum principle: 

To show that ut + f{x, t)ux = a{x, t)uxx satisfies a maximum principle we find a function 
that satisfies: 


F"{x) = F{x)/{2h{x)) 

where h{x) > a{x,t), and b{x) > cx^ + 1. We then consider the function z{x,t) = 
e^{F{x) + d) and following |^] we prove the maximum principle. 

We make use of general results about the solutions of F”{x) = F{x)/{2h{x)) in the 
proof. Because b{x) > cx^ + l,c > 0, it is easy to show |Q, p. 106, question 35]that the 
solution of the ODE must look like a solution of F"{x) = 0 for large enough x; i.e. F{x) 
looks like r±x + It is easy to see that if we let E(xg) = 1 and F'{xq) = 0, then the 
solution of the ODE will be concave up near xg. As F''{x) is positive at xg, we find that 
to the left of xg, the function F{x) must be positive and its derivative must be negative. 
Similarly, to the right of xg F{x) must be positive and its derivative must be positive. 
Thus F{x) > 1 for all x. Thus, F''{x) is also positive everywhere. That implies that 
F'{x) increase from some negative value to some positive one. Since asymptotically F{x) 
is linear, we find that F'{x) must tend to some constant value as x tends to Too. In what 
follows, we shall point out E(x)’s dependence on xg by using the notation F^^ to denote 
the solution of our ODE with initial data given at xg. 

We note that if we define Gxq^ = Fxfix) + d, then for all positive d, Gxo,d satisfies the 
inequality G”^ ^(x) < Gx^^/{‘2b(x)). As Fl^.^ is bounded, as d increase j_/GxQ,d\ tends 
to zero uniformly in d. 

Consider the function z{x,t) = GxQ^d{x)e^ . We find that zt = z and Zxx = 

Thus, 


Zt + f{x,f)Zx = 


> 


Gxo,dix)(^^ + f{x,t)G'x^^fix)e^ 

{GxQ,d{x)/2 + 6(x)G"g^rf(x) + /(x, 

Cl(x, i)Zxx “1“ e* {GxQ,d{x)/2 + /(x, t)G^g^rf(x)) . 


If we make d large enough we can make G as much larger than G' as we please. Thus, 
we find that for sufficiently large d, zt + Zxf{x,t) > a{x,t)zxx- Moreover, 2;(x,0) tends to 
infinity linearly in x and exponentially in t. 

We note that we can use e^^ rather than e* by letting F solve the equation F"{x) = 
cE(x)/(25(x)). This is what allows us to state that E{t) may be exponential and need not 
be sub-exponential. 

Now we modify one of the standard proofs of the maximum principle for the heat 
equation |^, pp. 216-218]. We consider w{x,t) = u{x,t) — €z{x,t) where e > 0. Clearly, 
w satisfies wt -b f{x,t)wx < a{x,t)wxx- If we consider this equation on a finite interval, 
[xi,X 2 ], then w satisfies the maximum principle: 


w{x,t) <'niax{w{xi,t),w{x 2 ,t), max t(;(x,0)). 

xe(xi,x2) 
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As z{x,t) > 0, we find that max 2 .g(a,j 0) < sup 3 .«(x, 0 ). Also, as z{x,t) increase 

linearly in space and exponentially in time (with any desired exponent) and, by assump¬ 
tion, u{x,t) increase more slowly, it is clear that for any xq, for any x of sufficiently large 
magnitude, w{x,t) < sup 3 ,M(x, 0 ). In fact it is not necessary that u{x,t) be strictly sub- 
linear. It is sufficient that it be sub-linear on an infinite sequence of points for which ±oo 
are limit points. If |u(a:,f)| < D{x)e°‘^ for all t and an infinite sequence of values, {x*}, 
that has as limit points ±oo, then we will be able to find an infinite sequence of points, 
Xi, for which w{x,t) < sup,,, tt(x, 0 ). 

Putting all of this together, we find that for any e if the magnitude of xi and X 2 
are sufficiently large and xi and X 2 belong to the sequence of points on which «(x, t) is 
sublinear, then w{xi,t),w{x 2 ,t) < sup,^, u{x, 0). Finally we note that as e ^ 0, w{xQ,t) 
u{xo,t). Thus, we find that u{xo,t) < supu(x,0). As our argument in no way depends 
on xq, we find that sup,j, tt(x, f) < sup,j, tt(x, 0). A simple corollary of this is that if all 
the hypotheses above are met, then if ut + f{x,t)ux = c + a{x,t)uxx, then s{x,t) = 
u{x, t) — ct satisfies the maximum principle. Thus, u{x, t) < sup^, tt(x, 0) -|- ct. Finally, by 
also considering —u we find that inf^, «(x, 0 ) -|- ct < tt(x, t) < sup^, u(x, 0 ) -|- ct. 

Remark 6.1. We see that if it is known that \u{x, t)\ < D{x)e°‘^ for all t and for an infinite 
sequence of x’s that run to Too, then u{x,t) satisfies a maximum principle. If one knows 
that the growth of u{x, t) in time is only exponential, then in order for u to fail to satisfy 
a maximum principle, it is necessary that u{x, t) have growth in x that is faster than any 
sub-linear function. In particular, u{x,t) > Mx/ln{x) for all sufficiently large x. 


Remark 6.2. We show that having some condition on f{x,t) is necessary. Consider the 
equation ut = b{x)uxx- If we let v = Ux, then we find that vt = {h{x)vx)x We have already 
shown that if 6(x) > -|- 1, then the equation F"h{x) = F has a solution that grows 

linearly at Too and whose derivative, F'{x) increase from some value at —oo to some value 
at oo. Let H{x) = F'{x). We find that {b{x)H'{x)) = H{x). Thus, e^H{x) is a solution 
of the equation for v. This is a bounded solution of the equation that satisfies neither a 
maximum nor a minimum principle. Of course, we can rewrite the PDE for v in the form 
Vt — b'{x)vx = b{x)vxx- Thus we see that it is imperative that some conditions be placed 
on f(x,t). Clearly this “hyperbolic term” can destabilize the parabolic PDE. 

If one specializes to b{x) which are even, then it is easy to say more. For such b{x) it is 
easy to see that the solution of b{x)F" = F with initial data F{0) = 0, E'(0) = a, a / 0 is 
odd, linear at infinity, and F'{x) is even and tends to a nonzero constant, k, at Too. F'{x) 
will always be greater than or equal to a. If b{x) = cosh(x), then (using more results on 
the asymptotic behavior of ODEs) we also find that F{x) kx + I exponentially fast at 
Too. Consider the function e*F(x) — kx — 1. It is initially bounded (it even tends to 0 at 
Too), it is a solution of u{x,t) = b{x)uxx, and it does not satisfy a maximum principle. 
We note that a function that solves the heat equation cannot blow up in this fashion. 
The bounds on the solution and the fact that the solution was initially bounded would be 
enough to guarantee that the solution remained bounded. 
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